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Abstract 

A short summary of different approaches to the parton energy loss problem is given. A particular 
attention is paid to the differences between various models. A possible solution to the problem 
of distinguishing competing approaches is discussed. 



1. Introduction 

For more than 2 decades, one of the most important goals in high-energy nuclear physics 
has been the study of the quark-gluon plasma (QGP for short). The existence of this deconfined 
phase of quarks and gluons has been unequivocally shown in many Lattice QCD (Quantum 
Chromodynamics) studies (for a recent study see Ref. [IJ). In those studies, the phase transition 
(cross-over) temperature between the hadronic phase and the QGP phase is shown to be roughly 
about 200 MeV when 3 light species (m, d, s) of dynamic quarks are taken into account. This 
value of transition temperature is actually somewhat of a conundrum. On the one hand, this 
temperature, although exceeding two trillion Kelvin, is nonetheless low enough to be accessible 
in accelerator experiments. On the other hand, the energy scale corresponding to this temperature 
is too low for QCD to be perturbative, making analytic calculations difficult. 

Nevertheless, combining the results from lattice QCD calculations, insights from perturba- 
tive thermal QCD calculations and also thermodynamics, we do know much about qualitative 
features of this new phase of the matter. Accordingly, many researchers have proposed many 
different signals of the formation of the QGP in heavy ion collisions. Many of these proposed 
signals utilize the fact that the energy density of QGP is extremely high. AiT - 200 MeV, 
the energy density easily exceeds 2GeV/fm^. Systems created at the Relativistic Heavy Ion 
Collider (RHIC) can reach maximum temperature of about 500 MeV. At this temperature, the 
energy density of QGP (composed of gluons and u, d, s quarks) can be as high as 100 GeV/fm^. 

Not surprisingly, the three most prominent QGP signals that emerged from RHIC experi- 
ments - the strong elliptic flow, the quenching of the high-energy (jet-) particles and the emer- 
gence of the medium-generated photons - are all measures of the high energy density and corre- 
sponding high pressure. Among them, the jet-quenching phenomenon is perhaps the most direct 
observation of the high energy density. Intuitively, it is easy to see that it will take an extremely 
dense matter to stop a particle with an extremely high energy. In this proceeding, a short sum- 
mary of the various theoretical concepts that go into calculating the jet quenching (equivalently, 
parton energy loss) is presented. It is, of course, impossible to do justice to the vast amount of 
work performed by many different researchers in this short proceeding. What I will do mainly is 
to use the McGill-AMY approach [2, 3, 4, ^ 6^ i ,8, _9] that I am most familiar with as an illus- 
trative example and highlight the differences between this and other approaches where possible. 
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Interested readers are directed to a more comprehensive review already in print ifioll and let me 
apologize here to people whose interesting works I cannot fully cover here for the lack of space. 



2. Schematic Idea of Jet Quenching 





Figure 1 : Schematic views of hard processes in hadron-hadron collisions (a) and nucleus-nucleus collisions where QGP 
is expected to be created (b). 

The idea behind the jet-quenching phenomenon is rather intuitive. As depicted in Figure [I] 
partons with high pj (jets) are produced when two hard partons from colliding hadrons undergo 
a hard scattering. In hadron-hadron collisions (Figure [T]-a), the jets then propagate and evolve 
in vacuum until they produce showers of particles in the detector. On the other hand, in heavy 
ion collisions jets produced by relatively rare hard collisions must propagate and evolve within 
the hot and dense medium (QGP) that is created by the rest of the system (Figure [T]-b). By 
measuring the difference between the outcome of these two experiments, one can then learn 
about the properties of the medium. The simplest consequence of having a dense medium is the 
energy-loss of the propagating parton, that is, "jet-quenching". 

The usual way of showing the effect of energy-loss is to show the ratio of the momentum 
spectrum from pp scatterings and suitably normalized momentum spectrum from AA scatterings. 
This ratio is usually referred to as the "Raa'' and defined as 



Nv,n,ry(dN''/dpT) 



(1) 



where h denotes the observed hadron species and A^binary is the average number of hard binary 
collisions at the specified impact parameter. (As far as the present author can find out, the first use 
of Raa appeared in 1990 in publications by M. Gyulassy and M. Pliimer flT] and M. Gyulassy 
and X.-N. Wang [12].) Experimental data from PHENIX and STAR collaborations spectacularly 
confirm that jet quenching is real. There is a sharp reduction of high energy particles in the 
away-side of the trigger particle fl3^ and Raa has a surprisingly small value of about 0.2 for all 
measured range of the high pj in central collisions 11411 . 

Theoretically, the calculations of the pp spectrum and the AA spectrum proceed as follows. 
For the pp case, the jet cross-section is given by the following schematic perturbative-QCD 
formula 
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where faip is the parton distribution function, dcTab^cdldt is the parton level cross-section and 
D/,/c is the fragmentation function. For AA collisions, we have 



^= r r faiA{xa,Qf)fmi^b,Qf)^^^nxc^x',\T,unDhiACQ) 0) 

.^geometry Jabcd 

where faiA is the parton distribution function of the nucleus A and the conditional probability 
Vixc — > x'J^T, u'') contains the effect of the medium that changes the parton from c to c'. Here 
T and m'^ denote the temperature and the flow velocity profiles throughout the evolution. We 
also need to integrate over the nuclear and collision geometry. Since the geometry is more or 
less fixed once the impact parameter range is known, the main task of a theoretician is then to 
calculate !P(xc — > x'^\T, u'') given the profiles of the temperature and the flow velocity. 

There are, however, some caveats that go with Eq.©. The pp formula (|2| is firmly based on 
the factorization theorem in QCD. In contrast, the factorization theorem that would put Eq.© on 
a firm setting is not yet ful ly p roven. (Gelis, Lappi and Venugopalan have taken promising initial 
steps in this direction jTslllal-) 

Another caveat is that the medium is finite and it is evolving all the time. The lifetime of the 
QGP created in a heavy ion collision is about 5 fm/c. The size of the system is about 10 fm. The 
hydrodynamic time scale is about tq ^ 0.5 fm. The mean free path of a particle in a QGP is also 
of order 1 fm. None of these are very large or very small compared to the others. Full accounting 
of these similar yet different length and time scales is therefore not an easy task. Inevitably, one 
needs to make some assumptions and approximations. 

With these caveats, nearly every theoretical approach to the energy-loss assumes that the 
above formula ([Sjl is at least a good approximation. And that's what will be assumed in this 
paper as well. 



3. Different Approaches to Parton Energy Loss 

As mentioned above, the main task of a theoretician working on parton energy loss is to 
calculate the in-medium modification function 'P(xc — > x'JT,u^). Collisions with the thermal 
particles cause the changes in the propagating parton. Hence, the fundamental quantity to cal- 
culate is the scattering cross-sections and the associated collision rates. If perturbation theory 
is valid, this would be a relatively straightforward calculation at least at the leading order. In a 
hot and dense matter, this is no-longer true: The dispersion relationship of a particle is no longer 
that in the vacuum. The scattering potentials are screened. Extra divergences appear due to the 
frequent soft exchanges. All these conspire to make the loop expansion invalid. 

For elastic scatterings, one can still regard the tree diagrams as the leading order lilTiflsll . For 
radiational processes, the above complications make non-trivial resummations necessary even for 
the leading order One must include both the elastic and the inelastic scattering processes jsl 
for phenomenology. Here, for brevity, we only discuss the inelastic radiationalprocess. 

The first study of radiational energy loss was conducted in Refs \21]. This method 



is often referred to as the BDMPS-Z approach. The main thesis of this approach is as follows. 
Consider a medium where the temperature is high enough that perturbation theory is valid. In 
this case, soft exchanges between the medium and the propagating parton result in the radiation 
of a hard collinear gluon. At the same time, the effect of multiple collisions is reduced because 
within the coherence length (or formation time) of the emitted gluon, all soft scatterings basi- 
cally count as a single one. This Landau-Pomeranchuck-Migdal (LPM) effect then necessitates 
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the resummation of all diagrams depicted in in Figure|2]in calculating the leading order gluon ra- 
diation rate. In the original BDMPS-Z approach some simplifications and approximations were 




Figure 2: Schematic views of hard scattering processes in hadron-hadron collisions (a) and nucleus-nucleus collisions 
where QGP is expected to be created. 

made. Specifically, the medium was assumed to be composed of random static scatterers. Also, 
the thermal screening effect was only approximately taken into account using the approach ad- 
vocated in Ref.[22J. The analysis was then carried out in the deep LPM regime and multiple 



emission is treated with the Poisson ansatz II23I1 . 

Subsequent approaches to the energy loss calculations use more or less the same starting 
point. One still needs to calculate and sum the diagrams depicted in Figure |2] The difference 
between various approaches may be classified by the way each approach (i) treats the scattering 
centers, (ii) resums the diagrams, and (iii) deals with the evolving medium. These differences 
may be summarized as follows. The acronyms in the following list are the initials of the main au- 
thors except HT (Higher- Twist) and McGill-AMY. Also, the references given are not exhaustive 
but just a starting point for further reading. 

• Scattering centers 

- Heavy static scattering centers (Cold medium) - BDMPS HI], Zakharov Q, GLV, 



1124 12511 ASW 112611, 

- Dynamic scatterers (Hot medium) - McGill-AMY HH, DGLV SIH, WHDG 

ii 



- General nuclear medium with a short correlation length - HT j 3^ 31 ] 

• Resummation schemes 

- Sum over diagrams with all possible soft interactions - BDMPS, McGill-AMY 

- Path integral representation of hard parton propagation - ASW, Zakharov 

- Reaction operator method with opacity expansion - GLV, DGLV, WHDG 

• Evolution schemes 

- Poisson ansatz - BDMPS, ASW, GLV, DGLV, WHDG 

- Fokker-Planck equation - McGill-AMY 

- Modified DGLAP equation (momentum space evolution) - HT 
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In this proceeding, I am going to use the McGill-AMY approach as the illustrative example to 
show what is involved in calculating each of the above items, simply because that is the one I am 
most familiar with. 



4. McGill-AMY Approach 

When the temperature is high enough, the asymptotic freedom property of the QCD makes 
it possible to treat QGP within the perturbation theory. The usual loop expansion, however, is 
no longer valid. The coherence effect (the LPM effect) makes it necessary to resum an infi- 
nite number of generalized ladder diagrams shown in Figure [3]-a. Technically, this comes about 




Figure 3: Diagrammatic representation of the generalized ladder diagram resummation. 



because connecting any of the three hard lines with a soft space-like gluon line (labeled with 
gT) introduces a pair of pinching poles in the loop-frequency space as illustrated in Figure |3}b. 
The resulting pinching-pole singularity, when regulated by the hard-thermal-loop self-energies 
cancels the factors of the coupling constants introduced by the additional interaction vertices. 
Within thermal field theory, Arnold, Moore and Yaffe rigorously proved that resummation of 
these generalized ladder diagrams is ncessary to get the leading order result (see |2] and refer- 
ences therein). Luckily, the resummation of the leading order contribution organizes itself into 



Figure 4: Diagrammatic representation of the Schwinger-Dyson equation for the radiation vertices. Figure courtesy of 
G. Qin. 

the Schwinger-Dyson type equation for the radiation vertex. Figure |4] shows a diagrammatic 
representation of the SD equation; The corresponding integral equation is 

2h = i6E{h,p, k)F(h) + g^J ^ {('^^ - CA/2)[F(h) - F(h-^q J] 

+(CA/2)[F(h) - F(h+p q J] + (CA/2)[F(h) - F(h-(p-k) q J]) (4) 



where 
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(5) 



Here p is the original hard parton momentum and k is the momentum of the radiated hard gluon. 
The masses appearing in the above equation are the medium induced thermal masses. The 2-D 
vector h is defined to be h = (p x k) x ey where ey is the chosen longitudinal direction and is a 
measure of the acollinearity. The diff'erential cross-section 



2 

Debye 



q2(q2 + ml^^^^) 

is the result of the hard-thermal-loop resummation. The open end of the vertex is then closed off" 
with the bare vertex and the appropriate statistical factors to get the gluon radiation rate 

In this approach, the medium consists of fully dynamic thermal quarks and gluons. Further- 
more, at least within the perturbation theory, the radiation rate so calculated is fully leading order 
in thermal QCD which sums an infinite number of the generalized ladder diagrams. 

The original AMY calculation of the radiation rates was carried out in order to obtain the 
leading order transport coefficients in hot QGP within the kinetic theory. In Ref.pj this is gener- 
alized to the case of propagating hard parton and its kinetic theory equation. The initial momen- 
tum distributions of the hard partons now evolve according to the following set of Fokker-Planck 
equations 

^ V^'^P^'^^i^t P,,(p)^J^^2P,(p^k) ^^^^ , 

dt Jk '^^ dkdt ^ dkdt ^ dkdt dkdt 

(7) 



which also includes the effect of absorbing thermal energy from the medium. For the full phe- 
nomenological study, this has to be supplemented by the local temperature and the flow infor- 
mation by independent hydrodynamics calculations. Since there is nothing intrinsically boost- 
invariant about our formulation, there is no restriction on what the underlying soft matter evolu- 
tion should be. The rates appearing above then become time and space dependent through the 
local temperature and the flow velocity. Finally, the resulting parton distribution at the final time 
is convoluted with the geometry and the vacuum fragmentation function to yield the medium 
modified fragmentation function ^ [stl : 

b(z,Q)^ \d^s^-^^^^^ fdpf-{'P,,(pj;pdD,(z',Q) + 'Pg(pf;pdDg{z',Q)), (8) 

The resulting Rf^A for shown in Figure|5l-a which includes both the effects of the radiational 
energy loss and the collisional energy loss |5]. One thing to notice is that Us ~ 0.3 does not 
validate perturbative treatment since this implies g ^ 2. At this point, McGill-AMY approach 
becomes a phenomenological model. However, I believe that this is the best one can do with 
current analytical tools since it at least includes all dynamic effects such as the energy-momentum 
conservation, broadening, thermal push, flavor change, elastic and inelastic energy losses except 
the interference between the vacuum and the in-medium processes lilOll . The last one is currently 
being worked on by the members of the McGill-AMY team. 
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Figure 5 : The nuclear modification factor Raa for Tr*^ calculated within the McGill-AMY approach. The panel (a) is 
from Ref.| 5] and the panel (b) is a preliminary result from MARTINI - an event generator being developed based on 
the McGill-AMY energy loss mechanism. In panel (a), the red solid line is the full calculation. The green broken line 
includes the eft'ect of radiative energy loss only while the blue dash-dot line includes the eft'ect of the collisional energy 
loss only. In panel (b), the green line is the main result with both the colHsional and the radiational energy loss. Here 
we set Q-j = 0.27 and the underlying soft-evolution is taken from 3-t-l D hydrodynamics calculation of Nonaka and Bass 
(3^1 . The data points are taken from fl4ll . 



5. Discussions 

In the previous section, a particular approach - McGill-AMY - was described and shown to 
yield results that compare very favorably with the experimental data. As far as the nuclear mod- 
ification factor Raa is concerned, other approaches can equally well describe the experimental 
Raa with a reasonable set of assumptions 1 10]. In this sense, the theoretical efforts to understand 
the jet quenching has been successful. However, the question remains: What do we really learn 
about the medium when different approaches with different assumptions about the medium can 
describe the experimental data equally well? 

For instance, consider the temperature. Since we are dealing with a thermalized system, 
finding the value of the temperature goes a long way to characterize that system. It sounds like 
measuring temperature ought be a simple task. Even this, however, becomes a non-trivial task 
among different models. Temperature appears as an explicit parameter of the calculation only 
in some models. In others, the medium is treated as dense, but cold to simplify the calculation. 
The connection to the QGP property is then made through the transport coefficient q = //^//l 
where fi is the soft momentum exchange scale and A is the mean free path. When soft exchange 
dominates the dynamics of the propagating hard parton, this should not be a bad approximation. 
However, this is still an indirect connection since fi and A are independent free parameters instead 
of functions of and T. Combined with the differing treatments of the multiple emissions, it 
makes the determination of the temperature within such models rather fragile especially since 
the value of Raa seems to be only weakly dependent on the value of q llssll . 

Recently, a collective effort to resolve this and other issues has been initiated by the TEC- 
HQM collaboration |34]. Right now the collaboration is performing standardized tests to see 
where exactly the differences between models lie. Since the shape of Raa is so featureless, I 
believe that the true test of models will come later when other observables such as the hard photon 
spectrum are calculated coherently within each approaches (for instance the y spectrum liSi lssIl ). 
It is therefore quite encouraging that Monte-Carlo event generators based on the current works 

7 



on the energy-loss mechanisms start to appear on the scene. In these proceedings, some results 
from YaJEM ["3?], JEWEL fy^, and Q-PYTHIA (s?) are reported. In the remaining space, let 
me introduce the event generator being worked on by McGill-AMY team - MARTINI (Modular 
Algorithm for Relativistic Treatment of heavy loN Interactions). The first, very preliminary 
result from MARTINI is shown in Figure |5}b. MARTINI is a modular modification of PYTHIA 
8.1 to take into account the energy loss of hard partons before they can hadronize (fragment) into 
showers. Hopefully, many more tests of the McGill-AMY approach can be performed on many 
different hadronic and electro-magnetic probes with this new development. 
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